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ABSTRACT 

The temperature f (r) in a stellar corona is computed under the circumstances 
that energy is supplied outward from the base of the corona only by thermal conduction. 
The heat flow equation is solved analytically under a variety of circumstances. In a 
corona of very low density the energy consumed by expansion of the corona can be 
neglected and ~f' ( 64 T ” , as in Chapman's original static coronal 

model. The result is a supersonic stellar wind with a velocity V ( oo) of the 

same order as the gravitational escape velocity £ ^ w . In a corona with 

medium density and sufficiently low temperature that V ( oo) is small compared 

to W , a near region, in which T (r) c c f' extends for some 

distance outward from the star before the far region j T(r) oc r ~ 2/7 
takes over. The result is a supersonic stellar wind velocity V ( GCt') °f the same 
order as the characteristic thermal velocity C c at the base of the corona. 

In a corona which is exceedingly dense, an intermediate region in which T(r) oc r ~ 
appears between the near and the far regions, which has the result of extending to large 
distance the point at which the coronal expansion becomes supersonic. In a corona 

* This work was supported by the National Aeronautics and Space Administration under 
Grant N AS A-NsG -96-60. 



which is exceedingly hot ( C 0 ^ w J the expansion becomes so violent 

that thermal conduction becomes negligible and the behavior of the corona is approxi- 
mately adiabatic. 

It is shown that any effect which tends to reduce the thermal conduc- 
tivity of the coronal gases at large distance from the star has the effect of enhancing 
the velocity of the stellar wind. 

Comparison with Chamberlain's earlier discussion of the solution of 
the momentum and heat flow equations in his "solar breeze" model shows that he made 
two self-consistent errors in his assumption that the energy flux in the solar wind is 
identically zero and that the gas motion is adiabatic at large radial distances from 
the Sun. It is shown that neither assumption is correct in a corona of finite density. 

It is shown, however, that the analytical form ) " ( it) Q C 1 / suggested by 

Chamberlain is obtained in the limit as the density of the corona is made large 
without limit, in which case all motion in the corona approaches zero. 

Application of the solutions of the heat flow equation to the sun 
— assuming that the solar corona is heated solely by thermal conduction— show that 
at least under present conditions the solar corona and wind would lie in the middle 
ground between high and low density and temperature. Assuming that they have 
coronas heated solely by conduction it is suggested that some of the giant stars with 
the low gravitational escape velocities, may fall into the high density case, and 
certain dwarfs into the low density case. Some of the very active stars may fall into 
the high temperature quasi-adiabatic case. 



I. INTRODUCTION 


In a previous paper (Parker, 1963a hereafter referred to as Paper I) the 
mass and momentum conservation equations were solved for a stellar corona in which 
the temperature was taken to be a given function T(r) of radial distance from 
the star. The present paper goes on from there to consider the form of ~T" (*'0 

in a corona in which energy transport is limited to thermal conduction. The paper 
is written from a purely academic point of view for the purpose of exploring the presently 
unknown dynamical properties of a conductive stellar corona. It is an open question at 
the present time to what extent the outer solar corona is supplied by thermal conduction 
from the low corona, as opposed to direct heating of the outer corona by the dissipation 
of wave motion originating beneath the corona. 

In this connection it is easy to show that thermal conduction is 
probably important, but it has not yet been possible to establish whether wave dissipation 
is, or is not, important, too. To show that thermal conduction is at least important, 
note that the observed solar wind strength at the orbit of Earth (Shklovskii, 1960; 
Gringauz, et al, 1960; Bridge et al, 1962; Bonetti, et al, 1962; Neugebauer and Snyder, 
1962) is of the order of 500 km/sec and, say, 5 or more ions/cm . The solar corona 
is observed to be at least as active at low solar latitudes, whence comes the solar 
wind observed near the plane of the ecliptic, as it is at high latitudes. Hence an upper 
limit to the energy carried away by expansion of the corona into the solar wind is 
obtained if it is assumed that the solar wind has the observed strength everywhere around 
the sun that it is observed to have near the ecliptic. The result is an estimated efflux 
from the sun of 7 x 10 protons/sec in the solar wind. The energy (gravitational 

) _9 

consumed by each proton in the wind is 5 x 10 ergs. 
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so that the efflux of 7 x 10 protons/sec means an energy consumption of 3.5 x 10 

ergs/sec by the expanding corona. Radiation losses above the base of the corona are 

27 

somewhat less than this, so 3.5 x 10 ergs/sec is a rough order of magnitude figure 


for the energy that is supplied to the outer corona. Now thermal conduction is expected 

to supply ljCT.)(yT). across )T = & where the 

temperature is "To and the temperature gradient is (V I 3* • Suppose 

that & * 1 . Billings and Lilliequist (1963) suggest that (V"T")o 

may be of the order of 3° per km. Using (2) this leads to an energy flux by thermal 

conduction of 1.1 x 10^ ergs/sec if x 10^ °K and 6.2 x 10^ ergs/sec 

if *T = 2 x 10^ °K. Thus thermal conduction could be the sole source of 
' © 

energy to the expanding outer corona (see discussion in Parker, 1963b). But the 
question is clearly open. Contemporary theory of coronal heating by wave dissipation 
(see for instance Osterbrock, 1961; Whitaker, 1963) is not sufficiently quantitative 
to be able to help in the decision. The one remaining approach to the problem is 
to solve the momentum and heat flow equations with simple assumptions concerning 
the configuration and structure of the corona to see if the observed solar wind 


velocity and density at the orbit of Earth can be accounted for in a quantitative way 
by the coronal temperature and density observed at the sun. A rough numerical investi- 
gation was begun by de Jager (1962) using Chapman's temperature distribution 
(Chapman, 1957, 1959) for a static corona. Noble and Scarf have recently begun an 
investigation of the solar corona and solar wind by numerical methods, giving a 


proper simultaneous solution of the momentum and heat flow equations. Their first paper 
(Noble and Scarf, 1963) on the problem shows that no energy source other than thermal 
conduction seems to be required beyond a couple of solar radii, within the present 
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uncertainties in the observations. However they point out that further investigation 
into some of the evident complications should be carried out before any final conclusion 
is reached. 

Altogether, then, it is clear that thermal conduction is an important 
process in determining the dynamical behavior of the solar corona. It may be presumed, 
therefore, that thermal conduction plays an important role in the coronas of many k i nds 
of stars other than the sun. In consequence of the evident widespread importance of 
thermal conductivity in stellar coronas the present paper undertakes a general study 
of the dynamical properties of the hypothetical stellar corona in which thermal conduction 
is the only form of eneiigy transport beyond some given radial distance f" “ 3- 
The purpose will be to examine the various qualitative dynamical features of the con- 
duction corona under different circumstances of temperature and density. We shall 
be concerned more with limiting cases, to illustrate the various features, rather 
than extensive numerical results for any single model. Thus the present study will add 
little or no quantitative information to the question of the dominance of thermal con- 
duction in the solar corona, discussed above. Rather it is aimed at illustrating the 
various asymptotic classes of coronal behavior so that the position occupied by the 
solar corona may be seen in its proper perspective. The aim is also to illustrate the 
dynamical possibilities available to stars other than the sun, many of which must have 
coronas with rather different values of density and temperature. 

Consider a stellar coronq in which T (r) is assumed to be deter- 
mined by the stationary heat flow equation* 


* See discussion and derivation in Appendix I. 
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hih m £}- N - 


- Zhld*l (,) 

dr /V ®r j 


where N is the number of ions per cm and V is the velocity of expansion, 
determined from TO) by the methods outlined in I. Radiation losses are neglected 
(Chapman, 1957, 1959). It is assumed throughout that T(°o) = 0. There are 
obviously certain special cases one can imagine for which T(a>) may not 
be essentially zero, but they will have to be taken up elsewhere. The thermal 
conductivity is denoted by Kj ("7”) and the numerical coefficients on the 
right hand side of (1) are appropriate to fully ionized hydrogen, for which (Chapman, 
1954; Spitzer, 1956) 


lj(T) = 6 )C/0 T erj>s/crvv sec. °K. (2) 

As in the previous paper the discussion will be limited to coronas in which the thermal 
velocity is small compared to the gravitational escape velocity, so that the corona is 
tightly bound to the star by the gravitational field. 

Now the heat flow equation has been discussed previously by 
Chamberlain (1961) in connection with his own ideas of the dynamics of the solar 
corona. Unfortunately he based his numerical solutions on the assumption that the total 
energy ftey* from the corona to r •= ©o is exactly zero. That is to say, 


Chamberlain postulated that the expanding corona forms the perfect thermal insulator, 
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so that with 10^ °K in the corona and 0 °K in interstellar space there is no heat transfer 


between. He also postulated that the flow was exactly adiabatic at large f , 
rather than solving the heat flow and momentum equations there, so he did not discover 
the error of his basic assumption. The analytical solutions given in the present paper 
show that the energy flow to r s flo is not zero, with the result that the stellar 
corona heated by thermal conduction expands into space with supersonic velocity, 

i / i/a 

rather than with an evaporative velocity going to zero like !•/ O as 

suggested by Chamberlain. It is interesting to note that Chamberlain's 
velocity dependence turns up in the present paper as the limiting form of the supersonic 
solution of the momentum and heat flow equations when the ratio of density to thermal 
conductivity, n A , becomes large without limit. The way in which this limit 
comes about is as follows: For any finite N j ^ the expansion becomes super- 
sonic at some critical distance fl (see discussion in I) and remains supersonic 

everywhere beyond; the distance Q increases without limit as " P G&j 

so that the supersonic portion of the solution moves out of the picture, leaving behind 
the evaporative solution V oC 1/r^ discussed by Chamberlain. Of course, 
the velocity of expansion goes to zero at the same time that its form approaches 


1 / r 

"A 


Va 


, so it is probably just as meaningful to say that in the limit as 
Oo the corona approaches stasis. 


II. BASIC CONSIDERATION,^ THE SOLUTION OF THE CORONAL HEAT 

FLOW EQUATION 


In order to carry out an analytical solution of the heat flow equation 
(1) it is necessary to consider the limiting conditions under which its solution may be 
desired. The study of the momentum equation in Paper I for arbitrary T(r) 
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showed two opposite extremes for expansion of a stellar corona. The first extreme is 
that as the temperature at the base of a corona is increased, a point is reached at which 
the effective enthalpy of the gas becomes equal in magnitude to the negative gravita- 
tional energy. When this occurs the corona is no longer bound to the star by the 
gravitational field. The expansion starts with the velocity of sound at the base of the 
corona, and the equivalent de Laval nozzle (see Paper I) looses its throat. The important 
point is that the energy loss to the star becomes so great with increasing coronal 
temperature that the energy transported by thermal conduction becomes negligible. 

The heat flow eqn. (1) becomes more or less irrelevant. The expansion of the corona 
becomes approximately adiabatic, and this case has been dealt with at some length 
elsewhere (Parker, 1960, 1963b). Adiabatic coronal expansion has the general 
property, already described, of starting with supersonic velocity at the base of the 
corona and decelerating outward to r sk OO . Any enthalpy excess 
over the gravitational energy is converted to kinetic energy. This circumstance of 
an overheated corona may perhaps have some application to the coronal outburst 
from the sun following a solar flare. This has been dealt with elsewhere, too 
(Parker, 1961). Altogether, then, the limit of high coronal temperature does not 
involve the heat flow equation in any fundamental way and it has been dealt with 
previously. It will be discussed no further here. 

The second and opposite extreme is in the limit of low temperature, 
in which the corona becomes to tightly bound gravitationally that the energy con- 
sumed by its expansion is negligible compared to the energy transported by thermal 
conduction. In this case thermal conduction is the dominant factor in determining 


T(r) 


and the effects of the heat flow equation are most readily 
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illustrated. Thus it is with the lower coronal temperatures that we shall be most con- 
cerned here, and consequently it is advisable at this point to consider a little more 
closely how the energy O Cn.v. a* ( 6 M* M /a)"] consumed by 
expansion compares with the energy flow OfrttWE] by thermal 
duction. Let represent the mass of the star, a the radius of the star, 

and M the mass of a hydrogen atom. It may be seen from either eqn. (25) or 
from eqn. (42) of paper I that very roughly 


corona 


con- 


as 


2 n v. = Of- 

\ akT 

T, becomes small. In general 6 M* m/ 3 k "To is of the 


order of 10 or more, so that decreasing To to one half means a reduction 


in 




y 0 by a factor of 10 . On the other hand the thermal flux is proportional 

to T. , which means a reduction only by about a factor of 10. At lower 

temperatures V 0 , and hence Ne v ( GM# M / ^) decreases even 
more rapidly in comparison to ^ (T.) * T . Altogether, then, 

it is evident that in the limit of low coronal temperature T. , the energy 
consumed by coronal expansion may be made arbitrarily small compared to the energy 
flowing outward through the corona by the mechanism of thermal conduction. In the 
limit of small coronal temperature and density, the coronal temperature 
approaches the static temperature distribution 


T(r) = T (-1) 


V7 
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originally proposed by Chapman (1957, 1959) for a static model of the solar corona. 

To facilitate solution of the heat flow equation (1) it is convenient 
to express it in terms of the same reduced variables as employed in Paper I in the 
discussion of the momentum and mass conservation equations. Thus let 3 desig- 
nate the radial distance at which the heat flow equation is assumed to become valid 
and put £ = r/a .Let e 1 = 2 kT / M , w ^ — 6My-/a J 
r\ = N / N a . The subscript zero is used to denote the value at M ~ . 

Then write =• ^ ( "To ) and 

if(T) - *)> f(‘Vco X ) 

so that f ) — (c / C* ) for fully ionized hydrogen. The heat 

flow equation (1) may be written 



upon using the condition for mass conservation, eqn. (6) of Paper I, Solving the 
momentum equation (5) of Paper I for (c f \ ) d n f d , substituting 

into (3), and integrating yields the usual energy equation 



N.v.a 2 M_ (<v 

^.*T. \ 



- Ek 

If.aT 


( 4 ) 
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where F“ is the energy flow per steradian at s — O . The 
quantity a represents the order of magnitude of the energy flux 

transported by thermal conduction outward from s — 1 . It is assumed 

that the temperature falls to zero at § OO , i.e. C l (oo ) = O . 

III. LOW DENSITY AND LOW TEMPERATURE APPROXIMATIONS 
The simplest asymptotic class of solutions of the energy equation (4) 
is in the low density approximation. For any given value of T. , the 
density (\[ 0 may be made sufficiently small that any quantity multiplied by 

N. v. a x / ffo a maybe neglected. With f ' Co ) 

sr (jL 1 f C. X ) , appropriate for ionized hydrogen, and with c l (-0 
2 

= Co , integration of (4) yields 


.'(i) = 


V? 


(5) 


The energy flux is 


FL = 


ft 


(6) 


The expansion v (§>') of the corona follows now from eqns (39) and (40) of Paper I 
wi th fS= 2/7 . The velocities and v(co) follow from eqns. 

(42), (43) and (45), (46), and are illustrated in Figure 4 of Paper I. As already noted, 
this low density approximation was applied by de Jager (1962) by numerical methods 
to the expansion of the solar corona. It is based on the assumption that the expansion 
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energy is small compared to the conduction energy 


N, 


Vo a 


M [i + 


<r< FZ 


(7) 


which can be achieved for any T. by making No sufficiently smal I. 
Hence we refer to it as the low density approximation. 

The condition (7) can also be satisfied for any No by making To 
sufficiently small, but if this involves making "To so small that w z /c a 
>> 1 0 , then v l (oo) << Vi x and the much better low temperatu re 
approximation is available. With W / C sT >> 1 0 there is a region 
0- } 1 } near star w ^ ere 5 C. *( 2 and V 1 ( 

are both small compared to ^ . In this region the energy equation 

(4) may be approximated as 


/ c 1 \ /2 JL (sl') ¥ NoVca^M W * J_ _ 
(, aV Ji " Ij.ii; 1 


FI 


l^o & To • 


( 8 ) 


This equation is valid out to ^ ^ , where C 1 ( § ) becomes comparable 

to w ^ . We recall from eqn. (8) of Paper I that 

W X / I > and ✓ *(§) are all of the same order of magnitude at 
the critical point. Hence the distance § t to the critical point must be of the 
same general order of magnitude as ^ ^ 

Now the total energy flux FTa at § ss CO is 

made up of the convection of kinetic energy and of thermal conduction. Denoting 



-li- 


the latter by 


we 


have 


fZ = N» v z (cd) + F7 . 


(9) 


Let it be assumed that the temperature is sufficiently low that 


No*o a v 4- 


( 10 ) 


This is a weaker restriction on No and than (7). Since 

is presumably comparable to V *CS) at the critical point near , it follows 

that all the terms on the right hand side of (4) must be small compared to F~a a 
by the time is reached. Under these conditions (8) is valid for all £ > i 

The term in w x / ^ on the right hand side of (8) is not correct beyond f i } 

but it is negligible there so incorrectness causes no error. 

Integration of (8) subject to the boundary conditions that C l (4) * 
and £*(oo) — O yields 


■(f) = 


Q. , 1^2 

r s 


*/? 


(11) 


where 


Q = 7 No 


4l(.iT 


( 12 ) 



- 12 - 


We require that O — is d in order that Si} be real everywhere 

in (i- ( ) . The energy flux transported by thermal conduction to ^ — CO 

is readily shown to be 


fT r 2 aT, (I-Q) 


The conduction flux across i- l is readily shown to be 


F. = (i + 9) 


(13) 


(14) 


which is larger than by [ust the amount of energy No*® 3 M w 

consumed in lifting the gas out of the gravitational field of the star. In terms of the 
definition of Q the basic energy relation 


/Vo v. 3 x M W * = Q £® 3 ~To 


(15) 


serves to determine v a as a function of No and T. , or if Y, 
ii known from the momentum equation, it serves to determine Cp . Now as we 
have already noted, the low temperature approximation is based upon the inequality 
(10) which is a weaker condition at low coronal temperatures, VY / 1* i O t 

than (7). If the coronal temperature is not low, then of course Y v ( Qo) — O 0 
and the two are equivalent, requiring small N/ o . The condition (10) may be ex- 
pressed in terms of Q as 
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w 


v*(cd) 


» 


Q ^ 

Ir9 


(16) 


This inequality establishes an u pper limit on Q for the validity of the low 

temperature approximation. Note that only as W */(•* ' ^ can Cp — > jL , 

There are several implications of (10) and (16) that should be noted. 

The inequality (16) requires that Q ^ 1 . It is obvious from (15) that 

Q > O , so altogether Q < Q 4 1 . The non- 

vanishing of 6? follows from the fact that with To > o we must 


have 


nonvanishing. For if V, b O , then c l (S) =c S 


- V7 


and the momentum equations predict a nonvanishing y, , in contradiction to 
the assumption that 

Now it follows from (11) that the outward decline of C '{§) 
must have an effective exponent ^ somewhere between 2/7 and 4/7, with 

s 4/7 close to the star and ^ =s 2/7 far away. The relative importance 
of the two regions depends upon the value of Cp .It was shown in Paper I 
that neither V* nor / ( 00 ) is very sensitive to the value of in 

this range. Thus for a given T. , both V, and V(oo) are very 

approximately independent of Q . Thus for a given Q , (16) is 
principally a restriction on T. , placing an upper limit on T. for the 
validity of the low temperature approximation. Eqn. (15) then gives once 

T. has been chosen (so as to satisfy (16) of course). The maximum value of 

No occurs for the maximum £p permitted by (16) for the chosen I o . 
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The lower is X the closer may G} approach to 1.0. This may be stated 
a little differently by considering that as Cp " 1 • 1 , (16) requires that 

must become small. To satisfy (15) it is then necessary that hj 0 become large at 
a very rapid rate. For intermediate values of , say Cp = 0.5, 

making /V* large, requires through (15), that T. must become small, which 
automatically satisfies (16). Of course, a moderate No , requiring a moderate 
T through (15), may satisfy (16) just as well in this case. Altogether, the 

maximum value of N* for any particular temperature T. occurs when Cp 
is as large as (16) will permit. For moderate or low T> this means that Q 
is near 1.0 and C ( § ) declines like J for some distance out from 

the star. 

Consider now the relative magnitudes of the distance to 

the critical point, the distance to the point where C i (s; 

becomes comparable to w 1 / j£ , and the distance § t to the point 

where C 1 ( f ) , given by (11) flattens out from the ^ near 


the star to the % 


- 2/7 


at large radial distance. Then is given by 


l. (8) of Paper I, and and § x may be defined to be 
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!t was already noted that is comparable to Jg ^ because 

and w y $ are comparable at §f e , More precisely, if c 
then solution of eqn. (8) of Paper I yields 


£ m 
f, 


2 (Zifi) 


'/(<-# 

= 0 ( 1 ). 


To demonstrate the relative magnitude of § . , note from the definition of 


N. Qnl = o ( FT ) 


Since y*(oo") is of the same order as W ^ , we have 

N.».a' -jr H = Q ( N.v. a 1 Mjt ") 


3» / 


- 4 " 

3t 

Then if (10) is to be satisfied, it follows that .It 

follows from the definition of that if , then (l - Q") 

>> Q , which is evidently equivalent to the statement (16). Physically 

this means that the region (A , f 0 near the star in which 

Vfe'i c-Vt 

c «c § terminates well before the critical point. In the vicinity of the 


-Id- 


entical point and beyond, then 


( I) 2 


'1- Q 


*/7 


( 20 ) 


From eqn. (8) of Paper I it follows that 


ro ( 7* 2 ^ I 


JG c*y (/- 9 ) 


vs 


( 21 ) 


This distance becomes very large as Q -+1 . The expansion velocity at the 

critical point is, from eqn. (9) of Paper I 


v’do = cHt) 



Noting that V*£oo) is of the same order as v 
matter to show that (16) and 
is equivalent to 


1 ( § J) , it is a simple 

are equivalent, which in turn 




(23) 
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Consider now the expansion velocity v( ^) that results from the 
temperature distribution (11) of the low temperature limit. The first approximation, 
given by eqns. (13) and (22) of Paper I, will be sufficient for the present purposes. De- 
fine 1(1) 


as 


!(£) = 


J 


vj W 


(24) 


u 1 C z (u) 


with C *(i; given by (11). Then from eqn. (13) of Paper I, the expansion 
velocity in ^ ^ is 


. . 4/? 

,*(§)= V. 


1 + 


1 - Q 
9 


t' 

ex r[z i ct)l <25 > 


The function 1(f) is related to the incomplete beta function and is expressed 
in terms of hypergeometric functions in Appendix X. The density in 2 < S« 


is given as 


n (f) = ey[- 2 I ( f)]. 




Putting ^ in (24) and using (22) for v ( it follows that 


H ■ fey ^ 


(26) 
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wIth I ( f c ) given in Appendix I. This serves to determine the mass flow 

/Vo Vo 3 ^ and, upon substitution into (15), it relates Q to 

N o and "To . 

In the region beyond 5 e the expansion velocity follows from 
eqn. (22) of Paper I as 

V? 


4 *(S) = 3cHiJ - 2 (-j, + 

- ZjlI (l- L) + 4J(S) 


(27) 


where 


J(i) = f ^ 


(28) 


wt 


.h c'(f) given by (11). The function J(i) can be expressed in 
terms of hypergeometric functions and is given in Appendix I. Noting that (11) 
approximates to (20) everywhere in 5 ^ S e , it follows that 


vT(°°) 


- ¥ [ 


87 , * f I fc t** 


7/s 


(29) 


To illustrate the results of the formal calculations under the low density 

(Q= o) and low temperature ( O ^ ^ 4 ) approximations the 

c C v/ l) >/z 

position of the critical point and the velocities / C% * i i 

v(ij/c., *,(eo)/c. are computed from the temperature distribution 
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given by (11) and plotted in Fig. 1 as a function of Co* . The asymptotic 

density ^ is plotted in Fig. 2. The values Q = 0, 4/7, 1 .0 were 

chosen to illustrate the full range of variation of Cp . Since (23) must be 

satisfied, it is obvious that =1.0 has no physical significance. It is included 

only to exhibit the mathematical boundary of the domain of ^ t etc. Cp = 4/7 
is valid only for w y^Co 5 , and of course Q = 0 is valid for all 

W 1 / Co Z . For the intermediate case Q = 4/7, about half of the. 

thermal energy is consumed by coronal expansion, the other half flowing out to 
5 = CO , There is the interesting matter of the cross-over of the density curves 
in Fig. 2 in the vicinity of W ft* 5 or 6. At large vj £» 

the curve for Q - 1 lies below the curve for Q = 0 for the reason that 
cVS) declines outward more rapidly for Q = 1. At smaller 

w y t* X the curves are reversed for the two reasons that with Cp ~ A. 

(^ = 4 / 7 ) the minimum value of w C* for — j. (g 
as Z fS 


iven 


by eqn. (41) of Paper I) is approached more quickly and the density 
at the minimum value of W l / <L 0 (given as 1 ( 3 / { 4-/3 - Zfi, x J\ 

by eqn. (47) of Paper I) is larger. 

To illustrate the range of validity of the approximations used in this 
section for a star with one solar mass we have plotted the thermal energy flux 
Z 7 ergs/sec steradian in Fig. 3 in comparison with the 

energy |\/o a L Fl w 1 ergs/sec steradian consumed in lifting the expanding 
corona in the gravitational field of the star. The numerical values of V e 
(shown by the broken line in Fig. 3(a)) were computed for an isothermal corona 
(4 = o} because such values are typical for ^3 


anywhere between 
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zero and 4/7 (see Fig. 4, Paper I). Fig. 3(a) was computed for as 1 f? 
for which the maximum temperature is 5.62 x 10^ °K ^ and Fig. 3(b) for 
b. — 4 , for which the maximum temperature is 1.4 x 10^ °K. A typical 

density in the lower solar corona, r = If?® is 10^ atoms/cm^. It is readily 
seen from Fig. 3(a) that the energy consumed by expansion is larger than that transported 
by thermal conductivity for T. * x 10^ °K, i.e. (7) is satisfied only for "To 

less than about 0.7 x 10 °K, On the other hand, a star with a coronal density of 
10 /cm will fall within the low density approximation for all I o . 

If it is assumed that the corona is heated by thermal conduction only 

beyond 4 f? © , then Fig. 3(b) is applicable. The density of the solar corona 

5 3 

at 4ff© is observed to be of the order of 10 atoms/cm (Van de Hu 1st, 1953), 
which leads to the conclusion that the energy consumed by expansion may be neglected 
if To £ 0.3 x 10 6 °K. 

It is evident from all this that the energy consumed by expansion of 
the solar corona is of the same order rather than much smaller, than the thermal energy 
flux transported by conduction. The solar corona seems to be too dense to fall into 
the low density limit, and somewhat too hot to fall into the low temperature limit. 

On the other hand, we have already pointed out that the solar corona is not so hot 
that it falls into the high temperature adiabatic limit. In the next section we take 
up the high density approximation. 

IV. HIGH DENSITY APPROXIMATION 


In this section we undertake the integration of the heat flow equa- 
tion (4) under the circumstances that No is taken to be extremely large while 
T. is maintained at some moderate value. The effect of making No 
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large is to make v ( 1} small, so that the high density approximation consists of 
integrating (4) neglecting only the term ^ f ) on the right hand side. It con- 
stitutes the next order of analytical complexity after the low temperature approximation. 
As we shall see, the solar corona is not sufficiently dense as to allow application of 
the high density approximation, but other stars with somewhat denser coronas may well 
fall into the high density category.* The principal interest in the high density approxi- 
mation is that it illustrates the behavior of the stellar corona in which the heat supply 
for a given temperature T. becomes small. The expansion velocity then 
exhibits a maximum at some intermediate distance and declines for a time before 
becoming supersonic at large distance. It is shown that in the limit as No 

c - '/z 

— ■ '■ ■ ■" » OO the expansion goes over into the algebraic form V (, f ) = Vj_ 3 
proposed by Chamberlain (1961) for the expansion of the solar corona, though of course 
Vj ■> O at the same time. 


1. Semi-Quantitative Discussion 

As a consequence of the greater complexity of the high density approxi- 
mation, there are a number of qualitative and semi-quantitative points that are best 
cleared up before presenting the formal solution of the energy equation. They are 
taken up under the paragraph headings which follow. 

(a) The Nonvanishing of h ^ : 

A fundamental point to be kept in mind throughout the discussion of 
the high density approximation is that is non zero for all finite N a . 


*lt should be noted that radiation losses become a serious concern when the density 
becomes, say, a factor of ten greater than the solar corona. This complication will 
not be included in the present formal solution of (4). 
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The nonvanishing of Fe is readily established by noting that for any No > 
no matter how large, the energy equation (4) may be reduced to 


-I'Gri®- *-£ 


T 

by choosing To small enough.* The reduction follows from the rapid decline 
of , discussed earlier. The boundary condition 

C — O means in this case that r m > o . We assume that increasing 

T. can only increase the total energy flow in the corona to § = CO , 
Hence 

FI > o 


for all higher F» . QED. As was pointed out earlier, Chamberlain (1961) 
made the error of supposing that was identically zero for all finite No 

and T . He made the supporting assumption that the coronal expansion 
becomes exactly adiabatic at large if . The low density and low temperature 
approximations have already demonstrated that the flow at large § is not at all 
adiabatic » . The same will again be true in the high density 

approximation. 


* This is a purely formal argument on the mathematical properties of eqn. (4), so 
the fact (4) may not be physically correct for a corona of arbitrarily low temperature 


is irrelevant. 
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(b) The Near Region: 

It has been shown that the energy flux is nonzero for all 

finite |\| 9 and vanishes only in the limit as N<> — ■> CO . On the other 
hand, the energy flux transported across § » 1 by thermal conduction is 
of the order of t. , which is essentially independent of No . 

With large /Vo only the small portion Fq» of A ~7~a is 

transported to J — Oo by conduction and convection. Hence, it is evident 
that most of a. / 9 must go into lifting the gas in the gravitational field 

of the star, which is N.v, a 1 M ** * . Hence, in order of magnitude 

(30) 

It follows that for a fixed T. , v . varies inversely with Ns and 

becomes arbitrarily small as /\/ a becomes large. The expansion velocity .(f) 
for i> 1 is proportional to , so that the velocity v ( ^ ) is 

generally inversely proportional to f\j a for fixed T. . It follows for 
sufficiently large N« , that in the region near the star both the term 
and the term in V *( §3 may be dropped from the right hand side of (4). The 
problem may be further simplified with the assumption that the temperature of the 
corona is fixed at some moderate value, say W c » ^ 1 O . Then 

for some distance out from the base of the corona the term C X ( f 3 ma y he 

neglected compared to w and (4) may be integrated to give 
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‘(f) £ c . 1 [l - Q (i - j,' 


V? 


(31) 


Now in this expression for < l (§) it must be required that Q 
is neither much greater than 1.0 , nor much less than 1.0. For if, on the one 

hand Q were much greater than 1.0, we would have «‘<f) 
vanishing near the star, which would require an enormous heat sink at the point of 
vanishing. No such sink exists. If on the other hand Cp were much less than 
1.0, would quickly level off at ^ Co i 

and a nearly isothermal coronal would be the result. But an isothermal corona yields 
a V 0 which is much too large to satisfy (30) as No — ^ CO . Thus, we 
conclude that, as was shown for the low temperature approximation, we must have 
Q 1 as No becomes large. It follows that 


e»(S) = ~ 


Vi 


(32) 


in the near region. The thermal flux from the base of the corona is then approximately 
415 . »t:/ 7 i and in place of (29) we may now write 


No a 1 M w 1 = y ifo ^ "X . 


(33) 


This expression determines V, as a function of No and "To . 

With the expression (32) for <*(!) in the near region it is 
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readily shown that 
side of (4) at § , 


5c*(f)/2 

where 


becomes equal to 


«ys 


on the right hand 



(34) 


The near region is defined to be (j- ; i i. ) 
equation for conservation of mass 


In it the velocity is given by the 



o(f)f 7 


(35) 


where n (. f J is determined by the hydrostatic barometric equation (see eqn.(12). 
Paper I) to be 



It follows at once that the coronal density and expansion velocity at the outer end 
of the near region are 



(37) 
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and 


gtv 2 

5 Co 2 -. 


V 7 A 


( 38 ) 


For fixed W 0 V / / Co 1 V ( Si ) is proportional to v Q , and hence inversely 
proportional to No in the near region. 

(c) The Intermediate Region: 

The energy relation (33) declares that V» is inversely proportional 
to No as N o becomes large while I „ is fixed. We know from Paper I 
that V 0 can be made small for fixed To only if the outward temperature 
decline in the corona becomes as steep as 1 / E . Th is is the role played 

by the intermediate region, which lies beyond the outer boundary S = S £ 
of the near region. The intermediate region, in which C * ( S) OC 1 / S y 
is the distinguishing feature of the high density approximation. We will have more to say 
on this later. 

In the vicinity of § , and beyond, the term 5 c 1 

cannot be neglected on the right hand side of (4). We obtain 


(fi) NkT) = [ £ <7(§) - 

VW Jfk>7 pi t 2 IJ. 


It is obvious from this differential equation that C (M) must decrease like 
1 . ^ Jr in the intermediate region. For if C Z ( decreased less 

rapidly, the right hand side would become positive and C Z would then 
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increase 


outward from the star. If decreased more rapidly, it would 

a ( c c - 4/7 

become negligible again, leading to C k $ J 5 which would be 

contrary to the assumption that C ( f decreased more rapidly than w § 
Now if C ^ ) is proportional to 4./ § , then the left hand side of the 


- s/z 


equation decreases like 5 , which rapidly becomes negligible compared 

c “ * 

to either term on the right hand side, both of which decline only as fast as S 
Equating the right hand side to zero, then, yields 


c l (f) = 


(39) 


throughout the intermediate region. It is this i-/ § temperature dependence 
through the intermediate region that leads to the small value of V* required by 
(33) for large N» . The gas flow through the intermediate region is approxi- 
mately adiabatic, since thermal conduction represented by the left hand side of (3), 
rapidly becomes negligible with increasing ^ .If v Z ( f ) is included 
on the right hand side of (4), the result is the wellknown adiabatic flow relation 


v‘(§) *%■ cYS) = o 


(40) 


If (39) were exact, rather than only approximate, the hydrostatic 
barometric equation would give 
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in the intermediate region. The equation for conservation of mass would then give 


(i) = v(S t )fe 


V* 


(42) 


The values of n jJ)and V (,^l) are given by (37) and (38) respectively. 

As a matter of fact, C ‘(f) must be a little less than (39) because of the neglect 
of the left hand side of (4) and because of the neglect of v C § ) . Thus 

must decrease a little more rapidly than (41), and the velocity is 
accordingly somewhat greater than given by (42). This is discussed further in para- 
graph (d). 

The intermediate region terminates where the individual terms on the 
right hand side of (4) become comparable to f~oa . If this occurs in the 

vicinity of , then in order of magnitude is defined by 


No v. a 1 M w- 

^3 



But the numerator of the left hand side of this equation is approximately equal to the 
energy flux F~^ transported across the base of the corona. Hence 
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(43) 


(d) The Far Region: 

Beyond f j , where fl can no longer be neglected on the 
right hand side of (4), the terms which become negligible are C and 

W j § , both of which vanish as § OO , Then if FL 

is written in the form (9), the energy equation becomes 




jjU*. 


= FL - 


« r< f M [v l (e°) -I'Hf)] . 


Since ^o* > O } it follows that F"[ and V (od) cannot both be identically 

equal to zero. Let it be supposed that F7 j ^ O . It is an easy matter to show 
from this that V ( ao} O / for neglecting V x , the energy equation yields 


C 



l 

Co 


JJ1 ) 

Z<.aT I J 


V? 


rv 1 
Co 



V? 
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It was shown in Paper I that for any outward temperature decline slower than -i/ ^ 
the expansion velocity approaches a constant v(ao) >0 as f becomes 

t - V? t 

large. The temperature here declines only as fast as 5 . Hence v (oo) 

is nonvanishing. The only other possibility is now that r.=o , but since 

O r this automatically gives v (go) ^ O . QED. Since 

^ s O j it follows that the expansion becomes supersonic in the far region. 
There must exist a critical point at which ( § ) 


crosses over 


cKl) • The temperature, given by (39) in the intermediate region is so high 
that the total enthalpy plus the gravitational energy of the gas is 


£ c l (S) -4- = o . 

2 ^ 


(44) 


We know from previous discussion that when this is the case, the expansion goes 
supersonic very soon after *'(S) declines less rapidly than Vi . Thus, 
we expect that in general order of magnitude 


S. « ° ( ■O 


(45) 


and 




( 46 ) 
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Presumably, then 


(») = o|V(fi)] 


(47) 


°(f) 


(48) 


It is immediately evident that the coefficient 2/5 on the right hand side of (39) may 
not be entirely correct toward the outer side of the intermediate region. It is readily seen 
from (40) that c*(S) is probably somewhat less than given by (39), with the 
consequence, mentioned earlier, that f\ ( § ) tends to be less than given by (41). 

The velocity v 1 ( §) tends to be larger than given by (42), so that it reaches a 
value ot C*C!,)] at The important point is that the small 

deviation from strictly adiabatic flow is an essential qualitative feature of the coronal 
expansion. 

(e) Discussion: 

The arguments presented in the foregoing paragraphs show that in the 
limit of large No the velocity Vo varies inversely with Na , as 
given by (33). With W * > > S C^/ 2* most of the energy is transported 
outward from § B 1 by thermal conduction. The convection of enthalpy may 
be neglected. The velocity increases throughout the near region (l f 
to a maximum v ( J given by (38) at the transition from the near to the 

intermediate region. Throughout the intermediate region the velocity declines. The 
conduction flux decreases outward through the near region like 1 f ^ 


whereas 
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i / — ^/? 

the convection of enthalpy decreases only as fast as -*■/ 5 • The near region 

ends where thermal conduction falls below convection. In the intermediate region 
thermal conduction may be neglected and the gas motion is approximately (but not 
exactly) adiabatic. The expansion velocity declines outward through the intermediate 

T 

region almost as rapidly as given by (42). The intermediate region ends at 
where the convection of energy becomes comparable to F~bo . The energy 
flux FZ> is nonvanishing for all finite Nio and guarantees that the expan- 
sion velocity will be supersonic at infinity. The smallness of the thermal conduction 
flux at suggests that the expansion proceeds to f = with only the 

total enthalpy that it has at . Consequently we suggest that the critical 

point lies not too far beyond fj and that v z (oa) is of the order of mag- 
nitude* of W 

It is of some interest to consider the form of v 1 C as 

(\ /«, — ^ GO . Eqn. (33) requires then V» — O like i/ . 

This reduction can only be brought about by the length § 3 of the intermediate 
region becoming large without limit. Then the flux of kinetic energy at § — O is 



(49) 


* It should be understood that the term "order of magnitude" as used through this paper 
is meant in the formal analytical sense. It is not meant necessarily to imply 
approximate arithmetic equality. 
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The thermal conduction flux at ^ — CD cannot be more than the conduction flux 
, which is smaller than the above expression by 1 ^ ^ 


at 


F t * O 



(50) 


The energy flux at ^ ■? CO is predominantly kinetic energy, and this kinetic energy 


goes to zero as N< 


CO . The result of N, 


CO is then that 


f — » CO so that the supersonic region beyond 3 > moves out of the pic- 

* 3 

ture. The approximate conditions given by (42) and (43) become exact and extend all 
the way to ^ « CO . We recognize this as the adiabatic subsonic coronal 
expansion proposed by Chamberlain (1961) as applying for all finite hi o . We 

have shown here that it arises only in the limit as hie — 05 . We note also, 

in this limit, that v ( §) — ^ O for all S / indicating that the corona approaches 
stasis. 


It is of interest to see under what circumstances the intermediate region, 
which is characteristic of the high density approximation, may be expected to occur 
in the corona of a star. There are a number of requirements for the existence of the 
intermediate region. In the first place the corona must be tightly bound to the star, 
or else the whole corona will be moving outward with supersonic velocity. This 
requires that the enthalpy plus gravitational energy in the low corona must be rather 
less than zero, 





( 51 ) 
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lt follows from (34) that this is just the requirement for the existence of the near region, 
i.e. ^ "> 1 .In addition to this, the intermediate region is subsonic, 

requiring that 

v(fj) C C f 1 ). (52) 


The velocity v ( § ^ J follows from (42), C (, ) follows from (34), and 

is given by (39). Thus (52) becomes 


N. » 


4K* 

i 

*35 a k Co 




(53) 


for the existence of the intermediate region. 

To see where the solar corona falls, put £ = 7 x 10^ cm, 

_ 15 2 2 _24 

V/ = 2.19 x 10 cm /sec , and put M = 2.06 x 10 gm for one atom in 

ten being helium. Then fora moderate temperature Vi = 10 ( = 1.6 

x 10^ °K) we have =25.4, = 1 .48 x 10 7 cm/sec, c(f*)=0.5 5 

x 10 7 cm/sec. The requirement (53) for the existence of the intermediate region becomes 


No » lMo“/tm 3 # 


. 8 3 

Observation suggests that Pv 0 is somewhere in the vicinity of 10 /cm . Raising 

T relaxes the requirement on No but we cannot go very far toward 

increasing T. because of the requirement (51). For the marginal case, that 
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y/ 1 / Co 1 is as small as 5 ( "To = 3'^X/0 °)C ~) , reached only 
under transient conditions at the peak of solar activity, we have § ^ = 5, Co 
2.1 x 10 7 cm/sec, C ) = 1.32 x 10 7 cm/sec and 


No » OlxlD*/ Cm 1 


8 3 

It is observed that No exceeds this value, but whether it exceeds 0.3 x 10 /cm 

an 

by a large enough amount for a long enough period of time to produce^intermediate 
region, we cannot say. It would be extremely interesting to see the phenomenon of an 
intermediate region in the solar corona sometime, somewhere beyond a distance of 
a few solar radii, but the numbers do not make it look very hopeful. Certainly any 
temporary intermediate region in the solar corona would be of limited radial extent, 
i.e. only a little larger than § ^ . 

It is not difficult to imagine stars with weaker gravitational fields 
and larger radii, so that * is smaller for a given value of 

, in whicl^the inequality (53) is easily satisfied. Perhaps some 
giant stars qualify for this condition. Then an intermediate region with a velocity 
maximum at § ^ would result. In some extreme case it is conceivable that N 0 
is so large that v ( CO ) might be only a few km/sec with the adiabatic solutions 


(42) and (43) extending far into space. 
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2. Formal Discussion 

With the substitutions 






+ 


JZ 

No vo a 1 Mw a 


No v* a 1 Hw 1 


Vs 



(54) 


(55) 


U = i 





Vi 


N o a v M W ^ 



(56) 


the energy equation (4) may be written in the reduced form 


y 



X- Y - U . 


(57) 



I 
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Consider the integration of (57) neglecting the velocity U . As § 
it is readily seen that approaches v wh« 


CO 


lere 


X. * 




s 7 6 


IT# a /o 


,V3 


FI 


5 Co y N^ 0 a"Mw a 


(58) 


Then since c^Coo) =■ O ■> it follows that 'NCX CO ) = o . Since 
FI approaches zero as No becomes large, it is evident that X » — > O , 

Thus the limiting solution, as No => CO , passes through the origin 

Y (0) = 0 . It is readily shown by repeated iteration of (57), taking advantage of the 
smallness of the left hand side, that the limiting solution is 

Y(X) = X - [x- (x-x s N - f X Nj 

* [ 1 - f (x -X - 1 x’ /z > f X "Yx - x **) s/! ] ■ ■ 


= X [l-X’'-5X 3 - ^x"No(X 4 )] 


(59) 


in the neighborhood of the origin. 
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When No is large but finite, Y( X^) i >s t . 

has the form 


Yw = (?r(x-x:r 


(60) 


as X ^ X x 

An important point to note before discussing Y ( Xj at large X 
is that the entire family of solutions Y (X,) tend converge rapidly toward the 
limiting solution ( IV« * — CD ^ as X increases from X <* .To 
demonstrate this convergence for small X « , consider a solution Y(X^) which 

lies a very small distance £ below the limitng solution Y^(X^ at X ‘ 5 - Xc 


i.e. 


Y(x.) = X(x.) - * ■ 


Let £ < < X £ ( X £ ) .In the neighborhood of fX., UX.)] 

let the solution be represented by 


Y(X) = X (x ) - ^ (X) 


so that € = 


K(X.) 


. Make X £ so smail that (59) reduces to 
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x(x;^x 


Then (57) reduces to 




) 


upon neglecting all terms second order in h . The desired solution of this 

equation is 


KX) = x* 


•3/2 




5/i 


v - 


sxr 


■f [ Js S S//2 «lt, 



£ £ 



- 1 ■=■ 




(61) 


Since X. <"< 1 , it is obvious that MX) decreases extremely 

rapidly with increasing ^ , demonstrating the convergence of the solution toward 

the limiting solution Ni ( for /No — Oo . 

Now consider the solution of (57) as ^ becomes large. This asymptotic 
limit js of physical interest because for moderate or low coronal temperature, 

1 and ^ becomes large toward the base of the corona, £ ~ 4 
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The asymptotic form of V" C X) ^ or l ar 9 e X is readily obtained from 

(57) by noting that X K ^ X . The result is 




(62) 


4 O 


This asymptotic form is valid for all large N/ e , as a consequence of the 
convergence of the solutions toward the limiting solution. 

The limiting solution was computed from (57) by numerical methods and is shown in 
Fig. 4. The other solutions of (57), which do not pass through the origin, are indicated in 
Fig. 4 by the short arrows which represent the slope J X / d X a * var ’ ous points' in 
the ^ X t Y ) plane. The convergence toward the limiting solution is clearly 
evident. 

The formal asymptotic form (62) reduces directly to (32) upon neglecting f~oo 
in (54). This gives formal proof of (33) for large No and ^ ~^/ £ <» 

Note that (54) - (56) can be considerably reduced through application of (33). The near 
region, then, is described by (62) and extends from X / - ^7X4) } corres- 
ponding to i = S, as defined in (34), to the base of the corona 

X. = S, 0/4) 

When X becomes small compared to X ^ , the solution lies close 

to the expansion of the limiting form given by (59), The leading term of this expansion 
gives simply X X , which reduces immediately to (39). The solution follows 
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y = X '"'Hi X nears X ^ . The intermediate region, then, 

lies between X us _) and X ± 

As X begins to approach X » the solution veers away from the limiting 
solution and follows the asymptotic form (60). With the aid of (33) we write 


X- X. 


i + 


FT 


FZ I 


(63) 


at 


is 


where, it will be remembered, is the total energy flux 4 A T./7 

% b= 4 . Substituting this into (60) and neglecting terms second order in 

i/S , it is readily shown that . C 1 ( § ) oC E • This 

constitutes the far region. 

Now consider the effect of inc luding U in (57). The effect of 0 
is to decrease jr/JX slightly. The decrease in eJ ^ / «=/ X 

negligible in the near region and throughout much of the intermediate region. The decrease 
becomes significant only in the outer portions of the intermediate region and in the far 
region. The result is that the solution Y ( X ) which includes U starts 

out at large X essentially coincident with the solution neglecting U . 

As the far region is approached, it declines less steeply than the solutions which neglect 
U . At small X the kinetic energy approaches a constant value, 

which when subtracted from F” ao leaves only the conduction flux 
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Thus X to must be modified by replacing FI on the right hand side of 

(48) with F7 . It was noted earlier that FT < Pols , with the result 

that X cp > which X vanishes, is moved much closer to the origin. The 
solution of (57) including U is sketched in Fig. 5. Noting that U O 

as No OO , ft is evident in this limit that the solution including Cj 

converges to the solution neglecting U , and both converge to the solution 
(59) through the origin. 

V. THE ROLE OF THE VARIATION OF (T) 

The discussion has thus far been confined to the hypothetical case of a 

stellar corona in which energy is transported beyond the base of the corona only by 

&/ s 

thermal conduction for which lj(T) is proportional to T . There are 
obviously many other possible forms for the conductivity >j(T) in the cir- 
cumstances encountered in stellar coronas. For instance transverse and disordered 

magnetic fields may cause *(T) to diminish outward from the star more rapidly 

--5/S 

than I . Or if we admit the possibility of un-ionized gas in, say, the outer 

atmosphere of a red giant, then /j(r; may decl ine more slowly than ”"T* 

— i — '/ 1 

being proportional only to - I . It is not possible at the present time to state 
precisely just what effects in *f(T) should occur in the solar corona, etc., 
except that there are several effects suggesting that t)(T) falls below the value 
given by (2) beyond some distance of the order of 0.1 a.u. or more. (See section VI.) 
Thus at the present time the best approach would seem to be an inquiry into the genera 
effects of deviation of f(T) from (2). Then when we understand the consequences 
l((T) in a general way it will be possible to state the results 


of variation in 
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of specific variations in *(T) when they become known. 

Several effects turn up in the general inquiry, the principal of which 
is that reduction in •SCT) beyond distances of several stellar radii may have 
the effect of enhancing coronal expansion. In fact it will be shown that a steady out- 
ward expansion of a stellar corona follows from the hydrodynamic equations if, and 
only if, 15 (T) decl ines outward from the star at a suitably rapid rate. 

1. Effects of the Form of (T) 

The simplest variation in *f( T) is a variation in the proportionality 
constant . It is readily seen from (4) that ^ a appears only in the 

ratio W./T. , so a reduction in ^ is equivalent to an increase in 

N 0 , which was discussed at some length in the preceding section. 

The next simplest variations in 5(T) are either to introduce a 
sharp cutoff in 1CT) at some fixed radial distance from the star or to assume 
that the functional dependence of on T is something other than 

T Vt , The effects of a cutoff in if(T) are considered under the 
next paragraph heading. The present discussion centers on the functional importance 


Let it be assumed that 


•j(T) = ij.fe.) < 64 > 

where V’" is a numerical constant, with the values ss O for uniform 

conductivity, T rz ± 

for an un-ionized gas, and rmS/t for a fully 
ionized gas. To demonstrate the necessary conditions on for coronal expansion 
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consider the case that the temperature C* or the density IV* at the base 
of the corona is small enough that the terms involving |\|» Vo on the right hand 
side of (4) may be neglected. Then subject to the boundary conditions that C ( l) 


= U' 


and 


c 1 (co) 


, integration yields 



i/ci+y; 


(65) 


and 



% a ck /» 


( 66 ) 


The discussion of the hydrodynamic momentum equation in Paper I shows that a corona 
will expand to supersonic velocity if,- or, ' ~ - . the temperature decreases outward 

less rapidly than J- / ^ .A more rapid decline gives a static corona. Thus it 
is necessary and sufficient that ^ ^ O in order for the corona to expand 
into a supersonic stellar wind. The corona is static if the conductivity is independent 
of the temperature or if the conductivity increases with decreasing temperature. 

There are some circumstances when fCT) may be more a function 
of position than of the local temperature. Then putting H CT) - r.f 
we have 


v 
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and 


f~Z — (l + &) § o * ~77 

for ^ y " <1 .* The corona will expand to yield a supersonic stellar wind only 

if % K O , i.e. only if 1^ declines outward from the star. 

The physical explanation for the outward declining thermal conductivity 
to favor coronal expansion is simple and straightforward. The coronal expansion leading 
to the stellar wind occurs in the region between the energy source at the base of the 
corona and the energy sink at Jj ** OO where the temperature vanishes. 

If the conductivity declines at large radial distance, the effect is to sever the con- 
nection with the cold § S CO . The result is a relatively slow outward 

temperature decline, leading to rapid coronal expansion. On the other hand, if the 
conductivity is large at large radial distance, the coronal thermal energy drains 
rapidly to ^ — Oo . The result is a rapid outward temperature decline, which 

discourages coronal expansion. Of course, the heat flow from the base of the corona 
is larger when the conductivity is large at large radial distance, but the flow passes 
to ^ — OO and does not serve to enhance coronal expansion. 

It would be a straightforward matter to repeat the calculations of the 
preceding sections for values of other than the ^ / “Z already 

considered, but the complete calculation is probably not of sufficient interest at the 

* If & K, ~ dL ^ the temperature falls to zero at finite § , at which 

point a heat sink is implied. This is physically unrealistic for the present discussion. 
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present time to justify presenting it here in its entirety. The iow density approximation 
has already been given. For the low temperature and high density approximations it 
is interesting to note that 


C 


z 





( 67 ) 


in the near region. The efflux of coronal gas is given by 


No 


Vo a 






t / c\ 4/(it y) 

in analogy to ( 33 ). In the far region C ( 3 J is proportional to 3 

which is just half as fast a decline as in the near region. Of course c * ( § ) 

is still proportional to VS in the intermediate region of the high density 

approximation. It is interesting to note that when ir 4 , the intermediate 

region disappears for the reason that C V ( declines at least as rapidly as 

wy s . Then with S £t» ' ^ Z W 1 at the base of the corona, 

5 c Ki)/Z remains less than w y § until the far region is reached where 

r. becomes non-negligible. 

The case ' *4 is of formal mathematical interest because the 

energy equation in the high density approximation (omitting V Z ( S) ) 


can be integrated in closed form (see Appendix II). When O the energy 
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equation can be integrated because V a is identically zero for the moderate and 

low temperatures 2^^ / 5 c, X >"> 1 considered here. 

2. Effects of a Cutoff of lf(T) 

There are a number of physical reasons, discussed a little later, why beyond 
some suitably large distance from the star ff ( T" ) may decline rather considerably 
below the value given by (2). If the decline should be abrupt, its effects can be 
approximated by introducing a sharp cutoff at some suitable radial distance 
up to which k ( ~T ) has its normal value and beyond which If (T ) is identically 

zero. Such an idealized form for (.T) lends itself to simple presentation of 

the physical consequences of a rapid decline in If (~7~) • 

It is evident from the previous discussion of the effects of the analytical 
form of Jf(T) that cutting off ty(T) at some distance may enhance 

coronal expansion. To illustrate the effects note that the motion of the coronal gas 
beyond s^. is completely adiabatic, because we are considering the hypothetical 
case that no energy source or sink besides thermal conduction is available above the 
base of the corona. The properties of an adiabatic corona have been discussed elsewhere 
(Parker, 1960). Briefly, the atmosphere will be static for moderate and low coronal 
temperatures unless 2 iv S C 0 : This may be shown by noting 


that the temperature in a static corona will be uniform out to 


; the adiabatic 


corona beyond can expand only if c x at the base of the adiabatic 

corona exceeds 2 w 

The really interesting case is when % ^ >> 2 tv 5~ C a 
in a corona of low density. Then in the low density approximation (dropping all terms 
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except FI on the right hand side of (4)) integration of eqn. (4) yields the 
temperature distribution 



and the energy flux 


f: = 4 <r.*T i - 

l i - i/f, 


7 A 


(70) 


at § » 1 . The energy is transported entirely by thermal conduction in the 

low corona and entirely by convection beyond , where C is 

obtained from the adiabatic condition V 1 
— V ( Oo) .At % - the energy flux has all been converted into 

kinetic energy. Thus 


FI = FT 4- M * z (co). 


(71) 


Now the velocity V ( oO ) can be computed from 
by the methods outlined in Paper I. The temperature c z ( will adjust 

itself so that v ( oc> ) in (71) will lead to an FI which is equal to the 

value given by (70). In the limit as No becomes small it is evident that F 7 
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must become small and C* ( ^ C* 1- , i.e. the corona becomes isothermal 

(1, !,) , with a temperature G* * . The dynamics of an isotherma 

adiabatic atmosphere have been treated elsewhere (Parker, 1960). A more interesting 
limit is obtained however if [\J * is fixed at some small value and § ^ is 
permitted to become large. Then C X ( § ) , as given by (69), declines rapidly from 
Co* to C*( close to the star and remains nearly isothermal with a 

temperature near c ( ) from there all the way out to . For an 

approximately isothermal atmosphere it Is readily shown from the momentum equation 
that 


V ( oo) ^ 



(72) 


for large 




We note that for a given value of 


velocity v ( co ) increases without bound as 


'5 


‘ <s,; 

becomes large. 


velocity , on the other hand is essentially independent of 3^ 

is large. Hence, it may be seen that with small /Vo 


r. - r. 


. Combining (70) and (71) yields 


, the 
The 
when 
(71) gives 


* The approximations involved are (a) that C ( $ ) is not precisely C (, 
near the star, where C ( %) — » C* and (b) C.(§3 is not precisely c(f^) 

as § becomes comparable to ^ because the convection term 
v V 5 c 1 " 2 ^ / § on the right hand side of (4) is not negligible 

there. Neither of these approximations is essential to the present argument, however. 
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No'/oa 1 J M v 1 (oo) 


* f Jf.aT. / 1 - 


c 4 (SJ 


W 


c.' 


(73) 


For any given C *” ( the left hand side of this equation can be made as 

large as desired by making 5^ sufficiently large. Thus, no matter how small 
may be M 0 , the left hand side of (73) can be made significantly greater than 

zero, with the result that on r '9^ band side is required by (73) 

to be significantly less than Co . Note however that C.*( VA l 

appears on the right hand side to such a high power that if C ( V is 

even ten percent less than Co , the term in C^( J ma y 

be neglected and 


/\/o vc a v \ !A v 2 (co) £ y ^ a~7l 


(74) 


This equation tells us that for small No , almost the entire energy flux 

^as § ^ bec omes I drq «. 

ZH.aTo/7 is converted into kinetic energy. Consequently, if 

f\j 0 is small, the resulting v ( oo ) may be enormous. Further, the 

velocity V * is determined principally by C 2 ( h) and diminishes 
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rapidly with declining 



, as noted earlier^ 



( 75 ) 


Then if 5^- is caused to become large without limit, V ( Os} in (72) 

becomes large without limit and C ( S ^ 3 diminishes in order to reduce y<» 

sufficiently to keep (74) satisfied. The point of this is that, within the framework of 
the formal heat flow equation (4), the cutoff of ^ (~T) at some extremely large 
radial distance can lead to an arbitrarily large stellar wind velocity 

V(oo) , This is true independently of the density No . The 

only effect of No is that for a larger N® , V 8 must be 

correspondingly smaller so that C ( § ^ ^ is smaller and ^ larger in 

order to achieve the same V ( 00 ) 


VI. SUMMARY AND DISCUSSION 


The present paper has concerned itself with the dynamics of the stellar 
corona in which heat is supplied above the base of the corona solely by thermal 
conduction. The corona of the sun may, or may not, be an example of pure conductive 
heating. One of the first points that should be taken up in the discussion of the 
formal mathematical examples presented above is the question of the general validity 
of eqn. (2) for the thermal conductivity. Eqn. (2) was derived for an infinitesimal 
heat flow in the direction parallel to magnetic field (Chapman and Cowling, 1958). 


Taking these two assumptions one at a time, we note that an absolute upper limit 
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to the heat flux is obtained by assuming that all of the electrons are moving in the 
direction of the heat flow with the rms thermal velocity. The resulting heat flux 
would then be k"T /2 'j ( 3 k. ”T/ rr\ ) ergs/cm^ sec. This maximum 

energy transport may be rather small if the coronal density is very low, and eqn.(2) 
is not valid if the actual flux begins to approach the maximum. The heat flux in 
a stellar corona is of the general order of ^ (l~) ^ ~T ( r ) ergs/sec stera- 
dian, at a radial distance f provided that the corona is dense enough to 

transport it. We must require then that 


'/i 

|(T)rT(r) « -f- Mr) ^ 


(76) 


If it be assumed that (T) is given by (2)^the inequality may be written in the 
numerical form 


N(.r)» 4X/0 5 T Vj y cm t 


(77) 


This inequality must be satisfied or else the effective value of the thermal conductivity 
will be depressed below eqn. (2); the effective conductivity ^ ( T) will be 
enough so that (76) is satisfied. Now the numerical values indicate that except in 
an extremely tenuous corona (77) will be satisfied near the star. On the other hand 
unless TCO drops off as fast as 4. / C at large f , ( 77) will 

not be satisfied as f — OO because N (r) is asymptotically 
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proportional to 1/r 1 , These circumstances may be illustrated by the numerical 

values for the solar corona. In the low corona where C = 10^ km and 

6 4 3 

T - 10 °K the inequality is N ^ 6x10 /cm , which is satisfied by a 

13 

large margin. On the other hand, at the orbit of Earth where f' = 1.5 x 10 cm 
and (/) seems to be of the general order of 10^ °K (Bonetti, et al, 1962; 

Neugebauer and Snyder, 1962) the inequality becomes N *>*> 3 /tm ^ 
which is just barely satisfied, if at all. Observations (see for instance Bonetti, et al, 
1962) suggest that N = 2 - 10/cm . A slightly higher temperature of 
3 x 10 5 °K yields N » 24/cm^, which is not satisfied. The conclusion is that, 
in the absence of any other effects, the low density of the solar wind must lead to a 
depression of the thermal conductivity below eqn. (2) at least at some distance 
beyond the orbit of Earth. The conductivity may often be depressed even inside the 
orbit of Earth. More complete and quantitative temperature and density measurements 
of the interplanetary plasma will have to be carried out before more can be said. 

Whatever may be the circumstances for the sun, it is evident that some 
caution must be exercised in application of the formal examples worked out in the 
text. The low density approximation and the cutoff with § ^ arbitrarily large 
are particularly suspect when it comes to actual models of existing coronas. As 
was stated in the beginning, the purpose of the formal models has been to illustrate 
the properties of the formal heat flow and momentum equations in a hypothetical 
stellar corona, concerning which there have been so many mistaken ideas. 

The second effect on the thermal conductivity is the well known 


channeling by magnetic fields. The reduction of the effective thermal conductivity 
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in the direction perpendicular to the local magnetic field is by the factor 

4* ^ (Chapman and Cowling, 1958) when the electrons are res- 

ponsible for most of the thermal conduction. Here cJ ca ' s the electron cyclotron 

frequency and L oe. is the deflection time for Coulomb collisions (Spitzer, 

2 

1956). When this reduction is by more than a factor of 10 it can be shown 
(Rosenbluth and Kaufmann, 1958; Vaughn-Wil liams and Haas, 1961) that the ions take 
over from the electrons and the reduction is not as large as the factor (l f CO 

&a tv) . 

In the low solar corona where N ^ I O ^ ^ the time 

between electron Coulomb collisions is of the order of 10 ^ sec, whereas a magnetic 
field of 1 gauss yields a cyclotron frequency of the order of 10^ radians/sec. The 
reduction of the thermal conductivity perpendicular to the field is clearly enormous. 

The same is true at the orbit of Earth and beyond. Altogether then, we may conclude 
that in stellar coronas and stellar winds the flow of energy by thermal conduction 
is channelled almost entirely along the magnetic fields. It follows at once that if 
the lines of force should become sinuous and generally non-radial, as they are 
observed to be in interplanetary space (McCracken, 1962; Smith et al, 1963), the 
effective path length will be increased and the cross section decreased, with a 
corresponding reduction in heat transport. It follows that if the field should become 
completely ^organized, as it appears to be at some distance beyond the orbit of 
Earth (Meyer, et al, 1956), the flow of heat might be cut off altogether. The result 
of reducing or cutting off thermal conduction in interplanetary space was shown to 
enhance the expansion of the corona. 

Now to summarize the results of the formal calculations. Assuming that 
the thermal conductivity is of the general form given by (2) it was shown that the 
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temperature in a corona of extremely low density declines outward like f 4 
At higher densities but for low temperature, a near region where the temperature 
declines like r appears at the base of the corona and extends for some 


distance outward before the decline goes over to T 


. If the density 


is high and the temperature is not low, an intermediate region where the temperature 
declines like f ” ^ appears between the near r and far r ' 

regions. It was shown, from the fact that a sufficiently low coronal temperature 
, reduces all these cases to the first, that the energy flow to f = CO i$ non- 
vanishing. Hence in all cases the coronal expansion becomes supersonic at large 'T 
to form a stellar wind. It was shown that the tendency of thermal conductivity to 
decline with decreasing temperature plays an essential role in bringing about 
supersonic coronal expansion. Coronal expansion is generally enhanced by any 
mechanism which tends to decrease the effective thermal conductivity at large dis- 
tances from the star. It can be asserted that so long as thermal conductivity is 
present the temperature declines outward from a star enough less rapidly than 1/ C 
that a supersonic stellar wind is the result. Only in the limit of large N.e Act.) 
is there a possibility of limiting coronal expansion to subsonic velocities. This 
limiting case may perhaps prove to be of interest jn the expansion of the coronas 
of some red giants where the low coronal temperature leads to very small ^(7*3 
Numerical estimates for the solar corona suggest that it is neither so 
tenuous, nor so cool, now so dense that it can be approximated by any of the three 
cases cited above. Nor is it so hot as to approach the adiabatic case. Rather 
the sun seems to lie solidly in the middle ground where much, but not all, of the 
energy transported by thermal conduction is consumed in the expansion. The 
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corona is hot enough that the ultimate expansion velocity is not only supersonic but 
it is comparable to the gravitational escape velocity from the base of the corona. 

Thus, simple models, such as the isothermal corona, or the isothermal-adiabatic corona 
(Parker, 1958, 1960) may be used to fit empirical data, but numerical methods, such 
as employed by Noble and Scarf (1963^ are required to deduce quantitative conduction 
models of the solar corona and solar wind from the mass, momentum, and energy 
conservation equations. 


Several of the qualitative conclusions resulting from the high and low 
density approximations may prove to be of interest for understanding the changes 
in the solar corona and solar wind over the 1 1 or 22-year cycle of solar activity. 
For instance, fora fixed coronal temperature "To the stellar wind velocity 


v( oo) 


declines to zero in the limit as f S/. 


OO . On the 


other hand y { ao} becomes comparable to the gravitational escape velocity 


2 ^ Vv in the limit as becomes small. Note also that for 

fixed i'he stellar wind velocity varies approximately as Tc . 

Now it is observed that both the temperature and the density in the low solar 
corona tend to decline during the years of minimum solar activity, particularly at 
high solar latitudes. One expects a simultaneous decline in the density of the 
solar wind and a decline in the solar wind flux }\Jo Vo d 1 . But it 


is not at all clear to what extent the solar wind velocity will decline, because 
the decline of |\| 0 tends to increase the velocity and the decline of T. 
tends to decrease the velocity. There is the additional possibility that with the 
declining density the thermal conductivity may cut off more completely at, say, 
oneor two a.u., which would further enhance v ( CO ) 


. It is not 
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Inconceivable that the net effect could be an increase in the solar wind velocity 
sometime during the years of low solar activity and a decrease of velocity sometime 
during the years of high solar activity. Such an effect might seem surprising, but 
at the present moment we cannot rule it out. Observations over the next 22-year 
solar cycle would settle the matter. The fundamental theoretical question concerns 
the relation between No and T. in a stellar corona. Presumably No 
and T. are determined by the wave dissipation which heats the corona, and 
while considerable progress has been made in this field (see for instance Osterbrock, 
1961 and Whitaker, 1963 and references therein) there is not yet any quantitative 
result of which we can be sure. 

When we come to consider the coronas of stars which are different from 
the sun, only the most general speculations can be made (see for instance Parker, 
1960, 1963b). Let it be assumed that there are stars in which the energy transport 
outward from the base of the corona is principally thermal conduction. From the 
general restriction that W / £ o ^ 5 it follows that C# must 
be very small for giant stars, and probably very large for dwarfs. The total heat 
flow is proportional to ( "To ) 3 "T# , so that if Co ex' 1 / 3 . 

we have that the heat flow is proportional to 3 . This may be very 

small for giants and large for dwarfs, suggesting that coronal expansion in giants 
may tend to fall into the high density category, and dwarfs into the low density 
and/or the low temperature category. It would be expected that some of the par- 
ticularly active stars might fall into the high temperature category with a violent 
and nearly adiabatic expansion. 
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APPENDIX I 

In order to evaluate the integral defined in eqn. (24) for the temperature 
cXS) given by (1 1), put 

\ s _v/ f 

Then wi-lh "2 



I(S) = 


W 


(i- <?) 


V? 


i(l-Q)/Q 


c . 1 


S/7 


io/7 


a -q)/q 


( i + *) 


i/7 


Note that by expanding the integrand in ascending powers of £ 
shown that 


f- 


it can be 


-ptr jy* (**• -y) - + 1 . _ X 

when Xj y < 1 , Here F" represents the hypergeometric function. 

When y exceeds 1.0, it is convenient to use the well known relation 
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2. — rAjr\ i rAi 1 1 - y) ss 0^2 Hil2 

* r(Min) 



r(ir\ l i*in j nf 4 i “ yj. 


A similar expression may be used when X exceeds 1,0. With fi\ 2 - "i/7 
and r\ * 5 / 7 /it follows that 




1 . 

7 > 







- i , £ 


1 . 

7 -> 



when (1-Q) S/Q <i .If (l-Q)/(?<t 

*> 1. / then 


1(8) 


7* v 

■iu rni 1 h i mm i 



1 2 . 

-> 7 



1-Q 
<P / 







12 . 

7 J 



-jfe raV7;r(/ 2 /7) 1 

r(V7D J 
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Finally, when ^ 


US) - 


_ 7w 




r / 5 2 J 2 

f T'ji 7i 


_£ 

J-c? 


_ J. F(£ > | • , - _ 2 _ 

- V7 ^T>T;T-> 


For purposes of numerical computation note that 


4 -,-*)• ^ 


44 


% - 


27 t , & % 

-x + -2- "X 3 

5» 343 


- + 6*< 
60025 6722 & 


x 1 + 


> * 4 


h( 7>^J 7i ■') ' 1 - ^ x + 4S_ x 2 _ J22_ 

42 


9*5 1 


445$ 


*> 


Since according to (16) we have (1- <?)f./c? >> i ^ it follows 
that with £ s Of /( L - Cp)f e 


I(SJ 


~ 7w J 


5c s (1-9) 


7/7 




V? 


f 1 + OC* 



Q < (i - 9 ) 


for 


and 
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V? 


- -L 

io 


- a/L£ 

5 


ro/7j> 


(^/ ^ [14 OCOj] 


for (l -Q>^) < £j> . 


To evaluate the integral defined In eqn.(24) for the temperature 
given by (11), 


JCI)- J £ t £T’ 


put 


W =■ » (i -Q) 

Q 


Then 


J(S) 


Ca 


*(4-a0 


V? , i/$< 


9 


V? 


7 


ft. 
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|t follows at once that 

J(f)= 7</(i- Q 1 )' 7 


*. . x . . -Q— 

7 ’ 7 J (1-Q)L 


1 r f l . i i 9 il 

[v, l 7. 7 - ^ - a-o)iJ 


rv 


?c. ! (1-Q) V7 
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APPENDIX II 


When the thermal conductivity 


* CO 


is proportional to the 


first power of the temperature, the energy equation may be written 


where 


and 


with 


y 21. = X - )r 

7 Jx 


y 


(Al) 


> 


- Q 


(A2) 


, £ 


r. 


(A3) 


> - 


^ Q 


(A4) 


9 


/v, 






o 3 I a 


(A5) 
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To integrate this differential equation let y •= V x. and t ® trt X . 
The result may be written 



(A6) 


where d is the integration constant and 

V! £ 4 >] ; 

Vz - 7 [(4* A')* - X ] 


(A7) 



The solutions through the origin are C -O and =■ /». x. 

X , The solution of physical interest is the one for which C l (co) 
This requires that y vanish at X *= X*, where 


in which case 



(A8) 


C 



V t 4 ^ 


(A9) 


For a dense corona, 

C<c L 


X 


» 4 < 1 as a consequence of the smallness of 
and the solution rapidly approaches the limiting solution 


Hence 
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y * Vx X as 7C 


increases. 



In the neighborhood of , the solution approximates to 


y = [ z (> - *«)] 


1/2 


(All) 


i. _ - 1/2 

which is to say that C oC 3 in the far region. 

The boundary condition C '’(l) “ Co serves to determine 
for a given W and 'OCoo . The quantity Xu 

can be determined only by simultaneous solution with the momentum approximation. 
In the low density approximation this is a simple matter, as was pointed out in the 
text. In the high density approximation it may <- be effect^for moderate or low 
coronal temperatures, for in that case , x 0- , which is nonvanishing for all 

finite , goes to zero in the limit as t\J 0 — Oq and (A 6 ) reduces to 


y * /z *■ 


(A 12) 
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In the I imit as No — 1 > oo we have k. “ Cp and y w Q 

where C » ) § — 1 . Then with \ sr Gy> (^5 It is readily 

shown with the aid of (A7) that 


_ i 

Q “ (1 - 5 e .’ ‘/ 2»') V2 , < ,3 > 


which reduces to 


No 3 X M * * ^ J. 


04 ) 


for moderate or low coronal temperatures. This is in agreement with (68) in the text 
since in the present case. Combining (A4) and (AT 3) there results 

the useful relation 


tl- fl‘* ( u i‘)"] 


in the limit of large No 

The solution of the energy equation for — J. has the interesting 
property that the temperature is proportional to J /f in both the near region 
and the intermediate region. Thus 



(A15) 


cX!) 





I 


all the way to the neighborhood of Otto where the far region begins. This may 
be seen from (67) or from (A 12). It is readily shown from eqn. (12) of Paper I that 
throughout this entire region 


n 




1 

1 


(A 16) 


and 


v(f) ^ V 0 S (A 17) 


The velocity becomes comparable to 



in the general vicinity of 




S 


(A 18) 
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where 


1 

»ycS- s/2 , 


and must be included in the energy equation from on and out to S "".OO. 

The critical point lies at some radial distance of, the same order as and 

the expansion velocity is supersonic beyond. 

The solutions of (A5), with C given by (A9) are plotted in Fig. 6 
for the moderate coronal temperature )\ = 0.25, for which = 1.133, 

Vi = 0.883. The limiting case of very high coronal density corresponds 

to the solution through the origin, x» =0 • The value of £ 0 1 

for this case is 11.3. Solutions for X» 0 (finite fs/o ) are plotted 
in Fig. 6 to show their rapid convergence toward the limiting solution (A12) with 
increasing x. 
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FIGURE CAPTIONS 

[U 

Fig. 1. A plot of | tj V L (Co} /c* > ( f c )/Va J and 

in the low temperature approximation as a function of w / Co 
for the three cases Cy* =0, 4/7, and 1.0. The case Cp = 0 represents 
diversion of only a vanishing portion of the eneigy transported by thermal conduction 
into lifting the expanding corona in the gravitational field of the star; Q = ^ 
represents diversion of all but a vanishing portion of the energy transported by 
thermal conduction; Op = 4/7 represents an intermediate case where about half 
the conduction energy is diverted into coronal expansion. The break in the curve 
of (/o / ' for Q = A/7 results from the use of the asymptotic 

expression for w V ^ 10 for the purpose of illustrating the order of 

the approximation involved. 

Fig. 2. A plot of the asymptotic density A Cf)f for large ^ in the low 

1 / 1 

temperature approximation as a function of ^ / Ca , 

Fig. 3. A plot of the energy flux 2 /fo & 71 / 7 from the base of the corona 

in the low density approximation as a function of temperature T 0 in compari- 
son with the energy v 0 3 ^ M ^ * for various values of [\) c 

consumed in lifting an isothermal expanding corona of the same temperature T. 
in the gravitational field of the star. The numerical values apply to the sun with 
(«> a = 1 R© - 7 x 10^ km and (b) as 4 . 

Fig. 4. A plot of the limiting solution as N* ->oo . The short- 

lines indicate the tangents to the family of solutions of (57) (neglecting U ) 


through the positions of the lines. 



-72- 


Fig. 5. The light lines represent a sketch of the solutions of (57) neglecting [J , as 
approximated by (59) and (61), The heavy line illustrates how the solution of (57) 
including U cuts across those solutions and reaches the X-axis at a small 
value of X. 


Fig. 6. A plot of y (*- ) 

= 0.25 and for 


from (A5) and (A9) for the intermediate value 
Xoo = 0, 0.01, 0.02, etc. 
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